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Abstract 

The previously developed Flux- Vector- Splitting (FVS) method was formulated here 
for the first-order version of Z4 formalism. Then, the characteristics of this method 
formulated for Z4 formalism were studied by numerically analyzing the evolution of 
two types of black holes (free and stuffed). Finally, these numerical results from the 
FVS method were compared with those from the Local Lax Friedrichs and Modified 
Local Lax Friedrichs methods to reveal the dependency of numerical solution by Z4 
formalism on the choice of numerical scheme. 
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I. Introduction 

Numerical relativity has been developed to solve Einstein equations. Various for- 
malisms derived from (3+1) ADM (Arnowit-Deser-Meisner) formalism [l], [2], [3] have 
been proposed to solve Einstein equation robustly. Among them, Z4 formalism ^ is 
constructed to have strong hyperbolicity, and as a result, techniques in CFD (Compu- 
tational Fluid Dynamics) can be used to solve Z4 formalism. Bona has also reported 
numerical results of black hole's evolutions by Z4 formlaism with numerical schemes, 
such as LLP (Local Lax Friedrichs) and MLLF (Modified-Local Lax Friedrichs) meth- 
ods [5]. In this current study, the FVS (Flux- Vector- Splitting) method recently pro- 
posed by Steger- Warming 6 was formulated here for the first-order version of Z4 
formalism. Then, numerical results from this formulated FVS method were then com- 
pared with results from the LLF and MLLF methods. The advantage of the FVS 
method over the LLF and MLLF methods is that the FVS method decomposes all 
characteristic waves of the convective matrix (the LLF and MLLF methods decom- 
pose only the maximum and minimum characteristic waves). Finally, the dependency 
of Z4 formalism on the choice of the numerical scheme was evaluated based on numer- 
ical analyses of three-dimensional evolution of two types of black holes (free [4j and 
stuffed 0]). 

II. Z4 formalism 

The first-order version of Z4 formalism is written as follows [7], 




Bk Diij — Bi Dkij, (3) 



BMi - Bi'A, 



ni Uk ^ -DI r>k , 



S{Q), 



(6) 



(4) 



(5) 



(1) 



(2) 



where X'^ijis defined with ordering parameter ^ as 




+ \S''J [Ar + A - (1 - ^« - '2Z^] , 



2Z,] 



(7) 
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with the foUowing definitions: 



D, = D,k*', (8) 



= D\^. (9) 

In eqs. (l)-(3), Ai,Bk^ and Dkij are defined as 

A, = d,lna, (10) 
Bk' = dkP\ (11) 
Dkij = ^dk-Jij. (12) 

In eqs. (l)-(5) and (10), a is called the lapse and /3* is the shift vector. 
In eq. (12), the intrinsic curvature is then given by 

^t■y^J + di 7y] = -2aA'zi, (13) 
ic curvature. 

In eqs. (4)-(6), source terms S{Kij), S{0) and S{Zi) are those expressed by Bona [7]. 
In eq. (1), Q is the rate of temporal variation of a as 

dta = -aQ, (14) 

and is expressed in harmonic coordinates as [7] 

Q = ~a'^dk\na + f(a)(trK --mQ). (15) 
a 

In eq. (15), / = is geodesic slicing, / = 1 is harmonic slicing, and / = 2/a is 
" 1 + log " slicing. 

In eq. (2), is the rate of the temporal variation of /3' as 

dtp' = -aQ\ (16) 
and is expressed in harmonic coordinates as [7] 

Q' = '^9kf3* ~ a-t'^' (a^7.fc ~dk\n^^dk In a) . (17) 
III. Numerical scheme 



In the numerical scheme to formulate the FSV method with Z4 formulism, we set 
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/3* = 0, because a value of is typical in the evolution of a single black hole 
Consequently, 31 independent variables needed to be solved: '^={Aa:, Ay, A^, D^xx, 

Dxxy, Dxxz, Dxyy, Dxyz, DxzZ, Dyxx, Dyxy , Dyxz, Dyyy, Dyyz, Dyzz, DzxX, Dzxy, 

Dzxz, D zyy , Dzyz, Dzzz, ^xx, ^ xy , ^xz, ^yy: ^yz, Kzz, ©5 Zx, Zy, Zz) are solved. 

Equations (1), (3)-(6) are rewritten as follows by using three 31 x 31 square matrices, 
namely, M, L and TV: 

dt"^ + dxaW^ + dyuL'i/ + dzuN'^/ = S (18) 

In eq. (18), aM can be rewritten as follows by using the matrix R, which is composed 
of the lines of eigenvectors of aM* and the diagonal matrix A, whose diagonal elements 
are eigenvalues of aM*: 

aM = RAR~\ (19) 

where R~^ is the inverse matrix of R. The formulation of R, R~^ and A is relatively 
straightforward (elements of M, R, R~^ and A in eq. (19) are shown in the Appendix). 
In contrast, the formulation of the eigenvectors of aL* and aA''* is complex. To simplify 
this formulation, we therefore rewrote eq. (18) as follows 

dt"2( + dxaW^ + iPy,xdyaLr + iPz.xdzaN'W = S, (20) 

where ^ (^y ^AxjAz ,Dyyy ,Dyxy ,Dyyz ,DyxX jDyxZ ,Dyzz ,Dxyy -jD x Xy -jDxyZ ,DxXX , 

Dxxz ,Dxzz ,D zyy ,D zxy ■jD zyz ,D zxx ,D zxz ,D zzz ^-^yy ,-^xy ,Kyz ,Kxx ,Kxz ,Kzz ^Zy ,Zx ,Zz^ 

and >^ (^2 5^1/5^3; ,D ZZZ ,D zyz ,D zXZ ,D zyy ,D zxy ,D ZXX ,Dyzz ,Dyyz ,DyxZ ^Dyyy , 

Dyxy ,Dyxx ,Dxzz ,Dxyz ,Dxxz ,Dxyy ,Dxxy ,Dxxx ,^zz ,^yz ,^xz ^-^yy ,-^xy ,^xx ,Zz ,Zy ,Zx^ 

In eq. (20), tpy,x and tpz,x are matrices whose elements are constants, either or 1, 
and defined as 



^ = i:y,Xy = ^Z,xW (21) 

In eq. (20), we introduce L and y, which are obtained by exchanging the superscript 
and subscript x for y in L and '?/ . As a, result, we consider dyaL'f in eq. (20) instead 
of dyuL'^ in eq. (18). Eigenvalues and eigenvectors of aL* are obtained by exchanging 
superscript x for y of 7*-* in each element of 7? and A shown in the Appendix. The 
same exchange of superscripts is done for R~^ to obtain the inverse matrix composed 
of the lines of eigenvectors of aL*. As with aN"^ , we exchange the superscript and 
subscript x for z in A*' and , thus yielding aNW. The eigenvalues and eigenvectors 
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of the matrix TV* are obtained by exchanging the superscript x for z of 7*^ in each 
element of 7? and A. The FVS method discretizes dxOtM'i/ in eq. (18) as [6] 

(9,«M<^)^^:fi±v2_:^, (22) 

= \ {=^.+1 +■^^- (G,+i - GO} , (23) 
J?, = aiMi<^„ (24) 
G, EEi?,iA.ii?ri'2<'„ (25) 

where i corresponds to the i*** numerical grid on the x-axis and /S.Xi is the grid size 
given by /S.Xi = \xi+\ — Xi\. Similar to the descretization of eq. (18) shown in eqs. 
(22)-(25), descretization is done for dyaLY and dzQ-NW in eq. (20). 
In contrast to the FVS method, the HLL (Harten-Lee-vanLeer) method [5] requires 
only eigenvalues without eigenvectors. The numerical flux in eq. (22) is defined 

by the HLL method as 

^ SR,i+l/2^i - SL,i+l/2-^i+l + SL,i+l/2SR,i+l/2 {^t+l - "^i) 

■J^i+1/2 — , (^O) 

SR,i+l/2 — SL,i+l/2 

where R = Right and L — Left, and SR^i^i/2 > and < in eq. (26) are 

defined as 

SR,i+i/2 = max (Ai+i/a) . (27) 
SL,i+i/2 = min (Ai+1/2) . (28) 

In Z4 formalism with f3i — 0, \sfi\ — |s||. Consequently, in the case of zero shift (i.e., 
I3i = 0), the HLL method is the same as the LLF(Local Lax-Friedrichs) method. 
In the modified Local Lax Friedrichs (MLLF) method proposed by [5], the numerical 
flux in eq. (26) is replaced as follows: 

^^+l/2 = + l ^ (Sfl,> + l'^« + l + SL,^'^^)} ■ (29) 

Both the LLF method and MLLF method use the same descretization procedure for 
dyoL'W and d^aN'^ in eq. (18). Comparison of eq. (23) with eqs. (26) and (29) 
reveals that the FVS method only decomposes all of the characteristic waves of the 
convective matrix aM, whereas the LLF and MLLF methods decompose only the 
maximum and minimum characteristic waves. Consequently, the highest accuracy 
is expected when all the waves are considered. The FVS, LLF and MLLF meth- 
ods considered in this study are first-order schemes. Higher order schemes can be 
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formulated by applying the MUSCL (Monotone Upstream-centered Schemes for Con- 
servation Laws) method [9] to all schemes. In this study, however, we analyzed only 
the spatial and temporal first-order schemes to focus on the choice of numerical scheme 
itself. 



IV. Initial conditions for the evolution of a single black hole 



In isotropic coordinates, the initial conditions for the evolution of a single black hole 
are as follows [1] 

^-,,1^=0 = 0, 7,,jt=o = *''<5»j, (30) 
a\t=o = l, /3"|t=o=0, (31) 
e|t=o = 0, Z,\t=o = Q- (32) 



In eq. (30), is defined as g] 



M\ M 



64 



X 2l -2 



P<f, (33) 



where p is the distance from the center of the black hole (i.e., the origin in isotropic 
coordinates), and M is the mass of the black hole. In eq. (15), the parameter / for 
the black hole's evolutions is set to "l-flog" slicing as 

f = I- (34) 
a 

In eq. (7), ^ for the evolutions of single black hole is set to 

i = 0. (35) 

In this study, we considered two types of black holes: a free black hole and a stuffed 
black hole. The difference between these two types is their energy density (r) distri- 
bution inside the Schwartzchild radius (p = M/2) as follows. 

T = (Free black Hole) (36) 
7^(/3<Tr) -^ = (p> M/2) (Stuffed black hole) (37) 



4M2 ^'^ 2 
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V. Numerical analysis of the evolution of a single black hole 

In this section, the evolution of the two types of single black holes is numerically ana- 
lyzed using the FVS, LLP and MLLF methods. The numerical grid is a 101 x 101 x 101 
unequally sized grid in Cartesian coordinates — 50M < X,Y,Z < 50M as shown in 
Fig. 1. The grid is denser near the origin (i.e., center of the black hole). First, we 
analyzed the evolution of a single free black hole using the same m values in eq. (15) 
reported by Bona, namely, m = and m = —3. 

Evolution of a free black hole for m = 

Figures 2A, B and C show the lapse profiles along the x-axis in increments of At/M — 
0.5 using the FVS, LLF and MLLF methods, respectively. Similar to the result ob- 
tained by Bona [1], the lapse did not stabilize, regardless of the method used. However, 
the time that the lapse became unstable was later for the FVS method than for either 
the LLF or MLLF methods, and the dynamics of this instability differed for the FVS 
method. In summary, the lapse by FVS method goes forward by the blow-up, although 
lapses by LLF and MLLF methods never go forward since the reverse wave begins to 
propagate into the backward. 

The validity of the code was verified by comparing our numerical results obtained 
using the FVS method with those obtained using the MMC (Marquina solver [10] 
with Monotonic Centered slope limiter) method by Bona 7 . Figure 3 shows the lapse 
profiles along the s-axis at t/M = 0,0.5, 1.5,2 and t/M = 14. In general, the lapse 
at < t/A/ < 2 obtained using the FVS method is similar to that obtained using 
the MMC method. The difference near the origin at t/M = 0.5, 1.0 between the FVS 
and MMC method is caused by the rough grid [7] used in the MMC method. The 
differences in results between two methods were significant at t/Ad — 14. These results 
reveal that the dynamics of the lapse depends on the choice of the numerical scheme. 
Figures 4 and 5 show profiles of and 0, respectively, along the x-axis in increments 
of At/M = 0.5. Zx is antisymmetric on both sides of the origin, whereas B = is 
symmetric. Figure 6 shows trK profiles in increments of At/M — 0.5. From eq. (15), 
the instability of the lapse can be explained by the growth of the domain trK < in 
Fig. 6. 

According to numerical results obtained using the MMC method reported by Bona 
[3], the lapse will always become unstable when m — —3 in eq. (15). We therefore 
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analyzed the evolution of a free black hole for m = — 3 in eq. (15). 
Evolution of a free black hole for m = —3 

Figure 7 shows lapse profiles along the x-axis in increments of At/M = 0.5 obtained 
using the FVS and MLLF methods. The lapse obtained using the FVS method be- 
came unstable when t/M exceeded 3.0, whereas the lapse obtained using the MLLF 
method remained stable. At t/M = 3.0, the lapse goes backward and overshoots unity 
remarkably. These dynamics of the instability of the lapse obtained using the FVS 
method is caused exclusively by the trK + 3Q term in eq. (15). 

Evolution of a stuffed black hole for m = 

Finally, we consider the evolution of a single stuffed black hole for m = in eq. 
(15). Figure 8 shows lapse profiles along the a:-axis in increments of At/M = 0.5. 
The lapse obtained using the FVS method became unstable when t/M exceeded 5.5, 
whereas the lapse obtained using the MLLF method still goes forward. The reverse 
wave going backwards emerged at t/AI — 5.0, which is earlier than the time of the 
emergence of the reverse wave in the evolution of the free black hole (see preceding sec- 
tion). Furthermore, the "serrated" profiles obtained using the MLLF method around 
x/M = 0.5, which is the Schwartzchild radius, were not observed in profiles obtained 
using the FVS method. 

In summary, the differences between the FVS method and the LLF and MLLF meth- 
ods are significant in the analysis of the evolution of free or stuffed black holes. In the 
evolution of a free black hole, the value of m in eq. (15) to avoid instability in the lapse 
obtained using the MLLF method causes the lapse obtained using the FVS method to 
become unstable earlier. Therefore, in the evolution of either a free or stuffed black 
hole, suitable values of / and m in eq. (15) are needed for the FVS method to avoid 
instability of the lapse. 

VI. Conclusions 

In this study, we formulated the FVS method for the first-order version of Z4 for- 
malism, and analyzed its characteristics by solving the evolution of free and stuffed 
black holes. The numerical results were then compared with those from the LLF and 
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MLLF methods, revealing that (a) the dynamics of the lapse in the evolution of free 
and stuffed black holes by the FVS method significantly differ from those by either the 
LLF or MLLF method and (b) a suitable variable rate of the lapse must be determined 
for the chosen numerical scheme such that the lapse does not become unstable. 
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APPENDIX: Elements of Matrix, M, R, R"^ and A In eq. (19) for ^ = In eq. (7) 



Non-zero elemnts of matrix Mij, (0 <i,j< 30) are 



Mo,2i = 

Mo,26 = 
M21,0 -- 
M21,ll 

M22,6 -- 
M22,14 
M23,2 = 
-M23,16 
M24,6 -- 
M25,ll 

M25,19 
M27,6 = 
Af27,10 
M27,14 
M27,18 
M27,29 
M28,26 
Af29,25 



= /7"", Mo.22 = 2/7"^ Mo.23 = 2/7^^ Mo,24 = /l"", Mo,25 = 2/7^^ 

= /7'^^Mo,27 = -/m,M3,21 = 1,M4,22 = 1,M5,23 = 1, M6,24 = 1, M7,25 = 1, M8,26 = 1, 
= 1,M21.6 = 7''^M21,7 = 27^^M21.8 = 7'',M21.10 = -7^^ 
= -7^'^,M21,16 = -7''^ Af21,17 = -7''^ Af21,28 = -2,M22,1 = ^, 
= -7'"^,M22,7 = -7'"'^,-^W22,10 = 7'"^,M22,11 = ^7^"'^, ^22,13 = ^J^" , 
= ^7^^,M22,16 = ^7"^^) -^22,18 = — ^7"^! Af22,19 = — ^7^^: -^22,29 = —1, 
= ^,M23,7 = -7"''', Af23,8 = -7'''', ^23,11 = ^7""^, A^23,13 = -^7^^,-^^23,14 = -^7"'': 
= ^7'^^,M23,17 = 7^^^,-^23,18 = -^7^^,-^23,19 = ■^7^"' , -^23,30 = -1, 

XX 



2' ' 2 

XX Tijr XX Ti/f XZ -AJ XZ -KjT XX 

= 7 ,iW24,10 = — 7 ,-M24,13=— 7 ,-^24,18= 7 , -™2B,7 = 7 , 
= -^7"'"^,M2B,13 = ^7'"'',-^\^25,14 = -^7''^-^^25,16 = - ^7^""', -^^25,18 = -^7""", 
= ^')'^\M26,S = 7'"'^, -^26,14 = 7'"'', -^26,17 = -7"""^ , -^26,19 = -7'"'', 

= -^^y^^y + ^^'^^yy^ m27,7 = -2'f''i"' + 2'f''i'\ M27,8 = -7"^^7"^^ + 7'^"^7^^ 

xy xy XX yy -k/t xy xz xx yz ji^ xz yy , xy yz 

= ^ ^ 7''*',M27,ii =7 7 - 7 7 ,-'^27,13 = -7 7 + 7 7 ^ 

yz I xy zz 71^ 

= -7 7"^ +7 7 ,-^^27,: 

xz yii xu iiz -h r xz yz xy zz 71 *■ xx 

= -y 7"" - 7 ''7- , M27,19 = 7 7 -7 7 ,A'27,28 = -7 , 

= -7'''', M27,30 = M28,22 = 7"^", ^-^28, 23 = 7^^"' , ^^28,24 = 7'"', 

= 27''^, M28,26 = 7'''', -^28,27 = -1, -^29,22 = -7'''", A/"29,24 = -7"^^, 

= — 7'^'',M30,23 = — 7'''^ > -^30,25 = —7^^^,-^30,26 = — 7^^^. 



Non-zero elemnts of matrix Rij, (0 <i,j< 30) are 



-27- -27-^ „ 

-ftO,0 — , ^,1 — , -^0,3 — ,^,4 — , 

' /yXX ' ' /yXX ' /yXX ' ' ^XX 

= — ,-Ro,9 = — , 



/(-^ -r in) 



Mo,'. 



f(-2 -r ni) 



(-1 + /)V7 



= ,-Ro,29 = -V77^,-R0,30 = VJT 



.y^^^y^^ ^^y^^^ ^'^^^^^ ^^y^^' 

Ri,i = 2, i^l,3 = + 7-, j?i,4 = + Ri,8 = -^-^ - 7", -Ri,9 = ^— ^ - 

' ' ' ^XX ' ' ' /yXX ' ' ' /yXX ' ' ' /yXX 

xy XZ ^xy xy ^xy xz ^xy xy 

R2,o = 2, i^2,3 = ^^-^ - 7"^ -R2,4 = - 7^^ -R2,8 = -^-^ + 7^^ -R2,9 = -^^-^ + 
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^3,16 
^3,19 
^3,20 
Rz,22 

^3,26 
J?3,27 
J?3,29 
-R4,15 
-^4,19 
-R4,21 
J?4,23 
-R4,2B 
Ri,27 
RiM 
-^5,22 

Rqa '■ 

-Re, 20 
Re, 23 
Re, 26 
Rt,3 '■ 

R7,9 '■ 

-Rs.s - 



^ -l + (7"^)V(-(7"^)'+7"V^) 



^l^y^^yy 

R3,18 _^xx^^xy'^2 _j_ ^^aia; j2yj/t/ ' 

/_ xz\2 w I j xy\2 zz 
— (l ) 1 + 17 ) 7 



/(7"^)^(m - 2) + 7-7^^ {1 + /(I - m)} ^ ^ 

^^J _ Y^^^jXX ^J7^'^^_(^^.xy^2 ^xx^yy"^ ' ^'^"^ ^^_„r.T _ ^ _j_ ^xx^yy'j 

7"^^V7^^(~(7^'')"^ -|- 72:3;^!/^^ ' "^^'^"^ _ya;3- (^ya;a')2 _|_ (^^,xx y2„^yy ^ 



yxy^yy 



-,R3,25 = 



-1 + 



— -yxx (^yxy'^2 _^ f^yxx'^2^yy '■ ^ 

(7"'')V(-(7"'')^ + 7""7^^) + (-1 + /(-I + m))/(-l + /) ^ 



(-),^^)2yra _ (y^!/)2^zz 27'"''(7^^7^^ - 'y'^y^y^) 



-K3,30 — 



"^^'^^ 7^'^y'7'"'^(— (7^^)^ -|- ryxx^yy'^ ' 



— 1,-R4 17 — 77 ' .„ r,-R4,18 — 7 TTT • 

' ^ixxn^xy\2 — ^xx^yy\ ' (^xy\2 — ^xx^yy ■ 



~,^y 



1^ 



{^r^xyy2 _ r^xx r^yy ^ ^'^^ ^r^xx{^_{^xyy2 _|_ r^xx^yy^ ' 

_ -y^^y^^) y^^((-y^!/)2 _|_ ^xx^yy-^ _ 2y^^y^!'y!'^ 
. 00 = ^ 



ryxx ^^xx ^^xy^2 _|_ ^xx^yy'j ' 4, 22 ^xx ^r^xx ^ ^,Yxy^2 _|_ ^xx^yy^ 



(7"")' 



_ ^xx^yy'^ ' ^'^^ ^^aiy^ _ ryxx^\ 
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7, -^4,26 = 



/ 

^ ^ ^ ^ ^ ''yxx f^^^xy'^2 _ ^xx^yy^ ' 

y»^!'(_(y»^^)2 ^ y'ca'y^^) _ya:2r^^X3/^2 _|_ ^xx^yy^ _|_ 2y^»^y^!'y2'^ 

ya:a;^y,xa;^_^y,a!y^2 _j_ ^xx^yy'j ' -^4>28 7'"^^7'^"^(— (7^^)^ _|_ yaiXyj/!/^ ' 

1 7'"'' 
1,^?B,17 = ,R5,21 = . , 

^XX r\jXX ^XX 



= , -^5,23 — —, -^5,27 — 



^,xy 



/ ? -^5.28 — " , 5 

/^XX r\jXX r\jXX /^XX 



■ 1-^6,9 = ■,Re,12 = 1, 

^xx ' ' ^xx ' ' ' 



ryX 



D 

(^^xy'^2 ^xx ^yy ^ ^'^^ _ (^jxy'^2 j_ yxx ^yy ^ 

'yxx (7^^^)^ — 7^^7^^ 2(7'^''7'^^ — 7"=^7!'^) 

1^6,21 = — , , , r75 T,Re,22 = 



^,yxy^2 _ jxxyyy ' -""O.^i ^^xx^^^xy^2 _ yxx^yy^ 
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— (^^xy^2 _|_ yxx^yy 



, -Re, 24 



, -Re, 27 



r 



7 

' "^^'^5 _^ya;y^2 _|_ ^xx^yy ' 



(7 



,xz-^2 _ ^xx^zz 



(^^xyy2 _|_ ^xx^yy ' 

,-R7,4 = 7T7Tr'-R7,l 



^ya;a!^_^ya;j/^2 _j_ ^xx^yy'j 



ri/^xy^xz ^xx_yz\ 

^ ^ -2(7 ^^7 -7 7" ) 

' ^'^'^ ^ySJa; ^y,a;i;^2 _ ^xx^yy'^ '■ 



27= 
7" 



27^ 



, -R7,ii = ;^,-R7 



2/yiCX 



2ya;a: ' 
1 



/7^ 



I -R7,28 



— — ,-R8,5 = 1,^?8,8 = -7,^8,21 = 



1 1 

7=, -R8,27 = 



12 



^9,13 = 
-Rl6,7 = 
-R21,17 

-R21,19 

-R21,20 

-R21,21 

^21,23 

i?21,25 

-R21,26 

-R21,27 
-R21,29 
-R22,17 

-R22,19 

-R22,21 

i?22,23 

J?22,2B 

i?22,27 

^23,17 

J?23,23 

-R24,17 

-R24,20 

il24,22 

-R24,2B 

-R24,28 
-R25,22 



■ l,i?10,12 = = l,i?12,10 = 1,^13,9 = 1,-Rl4,8 = 1, 

1,^?16,6 = 1,^?17,5 = 1,-Rl8,4 = l,i?19,3 = 1,^20,2 = 1, 
, ^?21,18 



_ h'^'f/i-m' + 7""7"") + (-1 + /(-I + m))/(-l + /) 



_ (7!:)V^-(7!!)!7 



,-R21 



_j_ (^^xx'^2^yy ' -""^i. ^xx^_^^xy'j2 _|_ ^xx^yy^ ' 

, -^21,24 



_|_ ^xx^yy^ ' 



.^^V xy\2 ^ XX ^ yv 

_|_ ^xx^yy^ ' 

(7"'^)V(-(7"^)' + 7""7^'^) + (-1 + /(-I + m))/(-l + /) 

,^21,28 ~ 



-^xx ^^xy''j2 _^ ^^xx^2^yy ' • ° ^xx ^_ (^yxy'j2 _^ ^xx^yy"^ ' 

1, J?21,30 = 1, 



l^xx i^_(^^xy^2 _j_ ^aJX^yj/J. 

, -R22,20 = 



R22Ai, 



^^XXyy^y 



_^^xj/^2 _|_ ryxx^yy ' -""^^i^u ^^xj/j2 _ r^xx^yy ' 
-y^y |_(-y^^)2 _|_ |(^^!/)2 _|_ ^xx^yy-^ _ 2^^a:^xy^B 



, -R22,26 = 



-^22,24 = 



(7a:3/)2 _ ^xx^yy ' 



fyxx'y^y 



(^■yxy'^2 _ yXXjyy ' 



7^ 



(7XHP _ ^xxjyy ■ 



, -R22,28 = 



7 



^xx j^^^xy^2 — yxx^yyy 



1 /y'!'^ ^^J/ 
^=,^?23,21 = , J?23,22 = , 

' ' /^/XX ' ' ^XX ' 



/7 



1 ^==2 ^""y 

, -R23,27 = —, ^23,28 = — , 



^ ryXX 



— {■y^y)^ + 7^^7^!' ' 
2^^w,y^^ _ 27'""' 7^^ 

— (^■yxy^2 _j_ ^xx^yy 

^yXX 

— (^^xy^2 _|_ ^xx^yy ' 
2^^a-y^z _ 2^xx^i;z 



^x 
, -R24,18 

R24.2 



l^xx^^y 



, -R24,: 



J?24 ,: 



(7x2/^2 _ ryxx^yy 

(^^^)2 _ ^xx^zz 
_^^xy"J2 _|_ ryxx^yy ' 
^^xx^xz 
_(^xj/^2 _|_ ^xx^j/H 
7^^ 

— (7Xt/)2 _|_ ^xx^y 



R24A9 = 



^^xy^2 — yyxx^y 



I -R24,: 



, -R24,27 



— (7x1/^ _|_ ^xx^yy 
(^^^)2 _ ^^^^^^ 
— (7X^)2 _|_ ^xx^yy ' 



— (7X^)2 ^ ^xx^yy ' 
1,^^25,28 = l,i?26,21 = 1, i?26,27 = 1, ^?27,20 = 1, ^?27,26 = 1, 
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R28,0 = "~;^'^2S,1 ~ ^^^^5^28, 19 = 1,-R28,2B = 1,-^29,1 = 1, -R29,18 = 1, 
^29, 24 = 1,-^30,0 = 1,^30,17 = 1,^^30,23 = 1- 



Non-zero elemnts of matrix Ri/, (0 < i,j < 30) are 

Ro,2 = 2' ^0,13 = 2 (j-^^ t''^) '-^0,14 = 2 ( 



•"0,18- 2 — + T j'-«o,i9-2(^ — -j;;— + ^ 

-Ki,i - o>«i.i3 - 7J ,«i,i4-;7 — r^+7 



2' ^'^■^ 2 y 7^"^ / 2 y / 2 y 7 

1 _ 1 / (7^"')^ zz\ p-1 _ ^ p-1 _ , p-1 _ , 

-"■1,19 — o I 7 I I -"■2,20 — -'-)-"3,19 — -'-)-t4,18 — -"-i 



■ 7 



2 y 7^ 

-Rb,17 = Ij^S.lS = li^7,15 = li^8,14 = 1; ^9,13 = 1; 

^10,12 = 15^11,11 = 1; ^12,10 — 15^13,9 — 

R-1 - ^ R-1 _ 1 R-1 _ T""^ 

-"■14,2 - -^::^'-"l4,5 - -L,-Kl4,7 - 

R-1 _ 7^^ p-1 _ 7^^ o-l _ 
-^14,8 — ~77--'^14,ll — ~ 7rT7 ■> ■'^14,13 



■XX ' 



XX ' 

1 

XX ' 



R-1 _ 7^^ R-i _ 7^^ p-i _ 7^^ 

-"■14,14 - ^-^y -«14,16 - -"14,17 ~ ~~ 

R-1 _ 7"^ R-1 _ 7^^ R-1 

-«i4,i8 - -^;:^:«i4,i9 - -^;:^.-Ki4,3o 

R-1 - 1 R-1 - 1 R-1 - 

-"15.1 - ~^:^'-"l5,4 - ^'-"15,6 - — , 

1 1 _ 7^^ R-1 _ 

-"is.'!- - -"15,10 - -;:^:-KiB,ii - 

o-l _ 7^ 



R-1 _ r R-1 _ _7 R-1 

-"■15,13 — ~ o a^a; ' -"15,14 — ~o xa: ' 15, 
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2,ya:a; ' io.i'i 27^^ J-0,1D 27^^ 

-1 _ 1^' -1 _ 7'^"' -1 _ 1 

-^15,18 = ^:^=-Rl5,19 = -^—^,Rl5,29 = —, 

R-1 i_ R-1 _ (-i + /)7"^ R-1 _ i-i + fyr 

-"16,0- J^xa! ' -"16,1 - /(7^^)2 '-"16,2- /(-ya;i)2 ' 

R-1 _i R-1 _ (7"")^(-2 + m) - + m) 

-"■16,3 - J-,«16,6 - (^yxxy ' 

R-1 _ 2(7"^7""(-2 + m) - 7""7^-(-l + m)) i _ (-y-)2(-2 + m) - 7"V"(-1 + m) 

-"16,7 - (,yxx)2 ' -"16,8 — (^xx)2 

R-1 _ (7"'')'(2 - m) + 7""7^^(-l + m) i _ j-yj-':{2 - m) + 7""7^^(-l + m) 

-"l6,10 - (^^xxy 5 -"16,11- (^jxx^2 ' 

(^^=^^^yy _Y''7m-2 + m) „-i _ (7" V— 7"''7"')(-2 + m) 

-"16,13 — / — N-i 



(^11)2 > --16,14 j-^j. 



R-1 _ 7-^7^- (2 - m) + 7-7^- (-1 + m) 1 _ (7"T(2 - + 7" + m) 

-"16,16- (^a:a.^2 1 -"16,17- (^a;i)2 

R-1 _ (7-V^-7-V-)(-2 + m) 1 _ (7-7^--7-V-)(-2 + m) 

-"16,18- (^xx)2 >-"16,19- (■^xx)2 ' 
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m 1 _ 7^^(-2 + m) p-i _7"^(-2 + m) 

~r7i -0-16,29 — TTTVJ '-"-16,30 ^ 



-"17,13 - 4 1^ 



+ 7" 



1 / _ Xll _ xz 



+ 7"^j ,j?r7'i8 



1 (^^^yf 



4 



I -^17,26 



2v^ 



?5 -"-17,30 — o 



'^1 -"-18,13 



1 

2' 

1 I 'i^^'-f'^ 



+ 7 ,-Ri8.i9-T 7^+7 



-1 

18,14 - ^ 



1 



~5 1 -"18,24 — 



^> -n-18,25 



2V7^ 



7 



^19,2 



^ p-1 

) ^iq 1 



47^ 



) -^19, 18 — 



R 



) -"-19,19 



4jx 



7 



2^/7^ 



2^7^ 



= 1-^19, 



19,24 



F) -Rl9,2 



1 -Rm 2 



2^7 
1 



2y'/ya:x ' 19,27 



J ^19, 



I -^20, 7 



^xy^xz _ ^xx^yz _^ 
J -^20,8 



19,28 — 2' 



2-\/7^^ \p^f^ 2-^7^^ 

_(^a:j/-j2 _|_ ^xx^j/y _ (^^a^^^^) -|- ^xz ^yy _ ^xy ^yz 
1-^20,11 = o /^S- )-R20,13 = 



2-/7^ 



7"^7''"-7"V 

, -TL. 



2v^ 



) -"-20,16 — 



2^7^^ 

_^xy^^xz _|_ ^xx^yz 
2y^p^ 



P_i _-(7^^)^+7"V 
I -"20,17 — " 



2v^ 



(7"V)+7"^''^ _ -(7-7'^-) +7-^- 1 _ 1 

7n=^ 1^20,19 — r, 5^20,27 — TT 



2\/7 

"5-^20,29 = 

^ -1 
"1-^21,14 



yXX 

s/Xy 



) -^90 .^n — 



2y^^' -^20,30 2^/^ 



D-1 _ V y 

'2^^^' 21,17- 2 ' 

7"" 0-1 _ 1 0-1 _ Vt^ 

2;^^, -"21,26 - 2 '-"22,7 - 2~' 



4 



^22,1 



7 



4V7^ 



1 Rii /t 



4^^xx ' 



4a/7^ 



, i?9 



4\/7^ 



) ^99 9 



1 1 _ 1 / (7^-)^ \ 1 _i ( t!!!!! , .A 
-4, -"23,13 - 4 — :^+7 j ,-"23,14 - 4 — :^+7 j 



1 fi^^^yf 



, -"23, 19 — 



]^ /^^xy^xz 



15 



7^^ „-i 



~^ 1 -"-23,25 ^ "^"y^ii^i -n-23,26 



-1 



1 p-1 _ 1 /'7"''7" 



1 / 7"''7 



^yz \ p-i _ 1 A7 

■7 I I -^24,14 



4 \ 7 



,xx 



■+7"M 



1 / i-y^' 



+ 7^^ ) , -^24^22 = 



7^^ -1 

: -^24, 25 



2a/7^ 



2V7 



7-« _ _ 

■"■25,2 — ^^„,^ I -"-25,13 — 



?) -^24. 



-24,29 — 9) 



>xy , ^,a:z _^^^^yy _|_ y^y^y^ 



^^xx 



r" 0-1 7 



7 -^25,18 — 



5-^25,19 — " 



^V^T^^ ' 25,23 2y^^ 



= 5 -^25. 24 — 



"^^'-"25,2 
1 

2' 



= ) -"-25,27 



1 

2v'7^' 



1 _ -(7"")2 +7-^7^ 

5 -"26,8 — ■ 



-(7 ) +7 7 p-i _ -7 7 + 7 7 p 

' -"26,7 - ' -"^e.s ~ 2^7 

_ ^aiai^yj/ j^y jXZ _ ^xx ^yz _ (^^^^j/y-j _|_ ^ajji^yz 

2^7^^^ '-^26,11 = 2^7^ 5-^26,13 = 2^7^^ ■ 



2^7^ 



1 -^9fi 1 



2^7^ 



R-i _ (7^^) 

) -tt-26,17 — 



2V7^ 



1-^26,19 



2\/7^ 



^ XZ ^ yz ^ xy ^ zz 

7 7^^ -7 ^^7 1 

) -"-26,27 



1 -^26,29 



2^7^ 



2^/7^ 

I ^26.30 = " 



2^7^ 



I -"27,14 — 75~^==S^i -"27,17 — ~ 



2 ''--i'.ii 

7^" p-i _ 1 
2^/;:px' -"27,26 - 2' 

"T; ) ^28,11 — 



1 ^28,13 — 



4^/7^ 



.'R28!l4 - -4-7^' 



^ I t>~l 
~~j I -"28,18 



?) -"29,1 



7 



D-1 _ _ 

4yO:^' -"28,19 — — g 

1.2 



7 D-i _ 1 

I -"28,25 — o5 



7 



?) -^29.2 



7 



•-29, 7 — 



27^^y77^ 

FTT7>rV7^ 



V7(-(7"'^)'+7""7"^) (-2 + m) ^ V7 (-7-^7"- + 7--7'^-_H_-2 + m) 

2(-l + /)7"^^a/7^ ' 
V7(-(7"')'+7""7"i(-2 + m) 
2(— 1 + /)7^^,/7^ ' 

V7((7"")^-7"V") (-2 + m) 

2(-l + /)7xx^/::pF '-"29,1 

V7(-7""7^'' + 7"V) (-2 + m) ^ V7(-7"^7"'+7"'^7^')(-2 + m) 

2(-l + /)7^^-y7^ ' 2(-l + /)7^^-y7^ ' 



y7(7"V-7"V^)(-2 + m) 

2(— 1 + f)jxx^^xx 
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R 



-29,16 
'29,18 



^29,21 



'29,25 
1 



R 

^29,28 

-^29, 30 
-^30,0 - 



R 



R 



R 



R 



R 



R 



R 



R 



R 



2(-l + /)7^^^7^ ' ^^'^^ 2(-l + /)7^^y7^ 



^ V7(7-7^'^-7-V-)(-2 + m) ^_i ^ y7(7"V^-7"V)(-2 + m) 

2[—l + fj-yxx^^xx ' 29,19 2(— 1 + j^'j-yxx^^xx 

_ 1 R-1 _ p-1 _ 7^^ o-i 

- 2 '-"29,22 - — ,-"29,23 - —,-"29,2 



27^ 



7- p-1 

7 



'29,27 



I -"■29,29 



,^^,«29,26 - 

y7(-2 + m) 1 
"2(-l + /)y^''^'" 

777"-M + m) 
2(-l + j)T'^^pT 

= , -"■30,1 



2-/m 



—27^^ + 2/7^^ ' 

y77"''(-2 + m) 
2(-l + /)7"'"'V7^ 



= , -^30,2 



'30,6 ■ 
'30,8 ■ 
^30,11 
'30,14 
'30,17 
•30,19 = 
'30,22 
'30I26 
-30,29 



2a/7\/7^'""^°'^ 2v77"''^\/7^' 2v77'^'^\/7^' 
^ V7((7"'')'-7"V) (-2 + m) ^ V7(7"^7^^-7""7^^)(-2 + m) 

2(— 1 + /)7^^>/7^ ' ^^'"^ (— 1 + /)7^^a/7^ ' 

^ V7((7-)^-7-7-) (-2 + m) , ^ V7(-(7"'^)^+7""7''")(-2 + m) 
2(-\ Jf- f)yxx ' 2(-l + /)7^'"V^7^ 

^ V7(-(7-^7-)+7-7^-)(-2 + m) ^_i ^ (7^-7^^ - 7- y-) (-2 + m) 
2(— 1 + /)7^^y'7^ ' ^"'"'^ 2(— 1 + /)7^^-^/7^ ' 

^(^^^^v^_^^!^^^^)(-2 + m) 1 ^ ^ + y--yy^) (_2 + m) 

2(-l + /)7^^y7^ ' 2(-l + /)7^»^.y7^ ' 

V7(-(7"^f +7"V^)(-2 + m) 
2(-l + 



-H-30,27 



2(-l + /)7"=^V7^ 

-^30,24 

2-/m 



1 V7(-7"^7^^+7"V^)(-2 + m) 
'■"■30.18 2(— 1 + /)7^^./7^ ' 

^30, 21 — 2 ' 



;:j^>-«30,23 - -"30,24 



^-1 ^ r: 

2/ya;ii: ' 30,25 ^xx ' 



27 



—27^^ + 2/7^ 



:> ^3 



^/7(-2 + 

2(-i + /)Vr 



y/^^i/(_2 + m) p_i V77""(-2 + m) 

F , -"30,30 



2(-l + /)7"^'^V^7^' 2(-l + /)7»'^V^ 

Diagonal elements of Kn (0 < i < 30) are 



Ao,0 = 0, Al,l = 0, A2,2 = 0, A3,3 = 0, A4,4 = 0, A5,B = 0, 
A6,6 = 0, A7,7 = 0, A8,8 = 0, A9,9 = 0, Alo,10 = 0, All, 11 = 0, 
Al2,12 = 0, Al3,l3 = 0, Al4,14 = 0, Al6,lB = 0, Al6,16 = 0, 

Ai7,i7 = — a-v^T^, Ai8,i8 = —a-\/7^, ■'^19,19 = —01-^/7'^"', A20, 20 = — a-v/7^,-'^2i,2i = — a^ 

A22,22 = — a-v/7^, A23,23 = a-s/T^, A24,24 = a-v/T^, A25,25 = a-v/7^, A26,26 = Ol\/^f^ , 



l^nai = 0^7^^, A28,28 = a^T^, A29,29 = — a V/7^: ^30,30 = a^/J^y^, 
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Figure Captions 

FIG.l 101 X 101 X 101 numerical grids for calculation of the evolution of free or stuffed 
black holes. 

FIG.2-(A) Lapse profiles along x-axis using the FVS method in the evolution of the 

free black hole for m = by At/M = 0.5. 

FIG.2-(B) Lapse profiles along x-axis using the LLF method in the evolution of the 
free black hole for m = by At/M = 0.5. 

FIG.2-(C) Lapse profiles along a;-axis using the MLLF method in the evolution of the 
free black hole for m = by At/M = 0.5. 

FIG. 3 Lapse profiles along a-axis using the FVS method and MMC method in the 
evolution of the free black hole for m = at < t/M < 2 (by At/M = 0.5) and 
t/M = 14. 

FIG. 4 Zx profiles along a;-axis using the FVS method in the evolution of the free black 
hole for m = by At/M = 0.5. 

FIG. 5 profiles along x-axis using the FVS method in the evolution of the free black 

hole for m = by At/M = 0.5. 

FIG. 6 trK profiles along x-axis using the FVS method in the evolution of the free 
black hole for rn = by At/M = 0.5. 

FIG. 7 Lapse profiles along x-axis using the FVS method and MLLF method in the 
evolution of the free black hole for m = —3 by At/M = 0.5. 

FIG. 8 Lapse profiles along x-axis using the FVS method and MLLF method in the 
evolution of the stuffed black hole for m = by At/M = 0.5. 
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FIG.l 101 X 101 X 101 numerical grids for calculation of the evolution of free 

stuffed black holes. 
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FIG.2-(A) Lapse profiles along a-axis using the FVS method in the evolution of the 
free black hole for m = by At/M = 0.5. 
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Lapse profiles along x axis byAt=0.5M (0-1 6M) 




X/M 



FIG.2-(B) Lapse profiles along x-axis using the LLF method in the evolution of the 
free black hole for m = by At/M = 0.5. 
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Lapse profiles along x axis byAt=0.5M (0-1 6M) 




X/M 



FIG.2-(C) Lapse profiles along a-axis using the MLLF method in the evolution of the 
free black hole for m = by At/M = 0.5. 
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Lapse profiles along x axis byAt=0.5M (0-2M) and at t=14M 
t=0.5M 




X/M 



FIG. 3 Lapse profiles along a-axis using the FVS method and MMC method in the 
evolution of the free black hole for m = at < t/M < 2 (by At/M = 0.5) and 

t/M = 14. 
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Zx profiles along x axis byAt=0.5M (0-19M) 




-4 -2 2 4 

X/M 



FIG. 4 Zx profiles along a-axis using the FVS method in the evolution of the free 
black hole for m = by At/M = 0.5. 
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© profiles along x axis byAt=0.5M (0-19M) 




X/M 



FIG. 5 G profiles along a-axis using the FVS method in the evolution of the free 
black hole for m = by At/M = 0.5. 
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trK profiles along x axis byAt=0.5M (0-19M) 




Free black hole 
m=0 

FVS method 



3 

X/M 



FIG. 6 trK profiles along a-axis using the FVS method in the evolution of the free 
black hole for m = by At/M = 0.5. 
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FIG. 7 Lapse profiles along x-axis using the FVS method and MLLF method in the 
evolution of the free black hole for m = — 3 by At/M = 0.5. 
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Lapse profiles along x axis byAt=0.5M (0-5. 5M) 




X/M 



FIG. 8 Lapse profiles along a-axis using the FVS method and MLLF method in the 
evolution of the stuffed black hole for m = by At/M = 0.5. 



